Ordered Loop Current States in Bilayer Graphene 



Lijun Zhu, Vivek Aji, Chandra M. Varma 
Department of Physics and Astronomy, University of California, Riverside, CA 92521, 



USA 



(N 

O 

(N 



S: 

a- 
& 

o 



00 

o 

O 



X: 



The zero gap non-interacting electronic structure of bilayer graphene is unstable to infinitesimal 
interactions. Indeed a gap is found in experiments but none of the states proposed so far is con- 
sistent with its properties. Here we show that a Magneto-Electric (ME) state is consistent with 
the experimental observations. This state breaks time-reversal symmetry through a pair of spon- 
taneously generated current loops in each layer, and has odd-parity with respect to the two layers. 
We also show that such a state wins over the Haldane-type state with spontaneous quantum Hall 
effect (AHE) which has even parity, for large enough interactions and low enough temperature. 
We also propose further experiments to confirm this conclusion and the possibility that at higher 
temperatures the Haldane-type state is realized. 



I. INTRODUCTION 

The non-interacting electronic state in bilayer 
graphene with AB (Bernal) stacking. See Fig. 1, has 
a pair of degenerate valence and conduction bands at 
two momentum points (0, ±47r/3-\/3) in the Brillouin 
zone. The energy varies quadratically with the momen- 
tum about these points, in contrast with the linear dis- 
persion in the single-layer graphene With hopping 
energy between the stacking carbon atoms in differ- 
ent layers, the other two pairs of conduction and valence 
bands have energy at ±i_L at the degeneracy points. Such 
a state with chemical potential at the charge neutrality 
point is unstable to infinitesimal electron-electron inter- 
actions. In weak-coupling approximation, in which the 
interaction energies are small compared to t^, one may 
look for instabilities restricting the Hamiltonian to the 
set of lowest energy conduction and valence bands. In 
such a reduced basis, a wide variety of states. Charge 
Density Wave, Valley Density Wave, Spin-Density wave. 
Anomalous Hall State (AHE) j3-[Tol| , have been proposed 
as possible ground states. However, the conductance ex- 
periments |lll4l5l |. done with high mobility samples and 
with both a top and a bottom gate to ensure that the 
chemical potential is at the charge neutrality point 
show an insulating state with characteristics which are 
not met by the proposed states. [l7| 

The two most important experimental findings can be 
summarized are that 

(1) The state at charge neutrality shows a gap in two- 
terminal conductance measurements with a conductance 
G smaller than the limit of measurement < IQ^^e^ /h for 



above which 

'1/2 



voltage V < EgQ, 

G^iV^ -El,r'l\ V>Ego, (1) 
with Ego ~ 2meV. 

(2) In small magnetic fields, the gap increases monoton- 
ically with a field, the conductance is consistent with a 
gap 



El^ + ul, u:, = eB/m*c (2) 
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(l/20)m, the free-electron mass. 



FIG. 1: (Color online) (a) Schematic of the bilayer lattice 
(top view). Bonds in layers 1 and 2 are indicated by solid 
and dashed lines, respectively. Here, A atoms in layer 1 (Al) 
and B atoms in layer 2 (B2) are stacked on top of each other, 
(b) is a unit cell for layer 1. Here 5^ and 5^ represent the 
bonds between the nearest (nn) and the next nearest (nnn) 
neighbors, respectively. 



The Charge-Density wave, and Valley-Density States 
are not consistent with a gap; the AFM state has a gap 
above which the conductance has the form of Eq. ([Ij , but 
it does not have field dependence of the gap of the form 
^ , (except in the unlikely case that conductance of only 
one layer is being measured) . The AHE state does have a 
a gap in the bulk and whose variation with field is also of 
the form But it also has a surface band with a quan- 
tized Hall conductance of ne^ /h with n depending on the 
degeneracy, which in this case is 4. This contributes to 
the two probe conductance measurements and therefore 
is not consistent with the observations either. 

An AHE state is possible if either inversion symmetry 
is preserved or if it is broken, the product of inversion and 
time-reversal is also broken [3, [l^ . In bilayer graphene 
with AB stacking, the stacked atoms in each layer are at 
points of inversion. See Fig. ([T]). In this situation, AHE 
is possible only if two conditions are met: (1) the direc- 
tion of the flux in the triangle formed by the stacked and 
the unstacked atoms in each layer is the same, (the mag- 
nitude is in general may be different). (2) the direction 
of flux in the triangle formed by the stacked (unstacked) 
atoms in layer 1 is the same as that of the unstacked 
(stacked) atoms in layer 2. The other three possibilities 
violate inversion with the product of inversion and time- 
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reversal preserved and therefore do not have an AHE. 
In the weak coupling approximation (V^„„„ ^ t±, where 
Vnnn is the next neighbor interaction), AHE is the only 
state accessible among these four. We show that when 
the coupling is larger (Ki„„ > 0{t±)), the ground state is 
the ME state, with a flux pattern depicted in Fig. ([3ji-f) 
and has properties consistent with the experiments. 



we confine our attention to symmetry breaking due to 
Vnnn and assume that V„„ only renormalizes the param- 
eters used in the calculation. The mean field analysis of 
the model starts with the decomposition of the diagonal 
in spin part of the density interactions between sites as 
was done to derive time-reversal breaking ME states in 
cuprates [l^l 



II. MODEL 

The non-interacting model has a nearest neighbor hop- 
ping i, a next nearest neighbor hopping ti and the hop- 
ping between the stacked carbon atoms t± = tAi,B2- We 
consider a nearest neighbor interaction parameter Vnn 
and a next nearest neighbor interaction parameter Vnnn, 
which controls the many-body effects. On-site interac- 
tions only renormalize the parameters for the states we 
consider and therefore will not appear. Other possible 
terms in the kinetic energy such as tA2,Bi ~ 0.3eV (the 
trigonal warping term) do not make any essential differ- 
ence,. 

The Hamiltonian is 

H = Hi+H2 + Hr2 + H,nt (3) 
Hi2 ^ a|j62i + h.c. 

i 

Hint = Vnn ^ + Kmn ^ ''^H'^^li^S^- 

I is the layer index taking the values (1,2), nu = (^nii^n 
where c is the fermion operator, au or hu and nu the 
charge density operator. {5i to 5z) are the three vectors 
connecting the nearest neighbor sites on each layer and 
(^1 to 5q) are the 6 vectors connecting the next nearest 
neighbor sites on each layer. 

Estimates for the parameters extracted either from 
band structure calculations or deduced from experimen- 
tal measurements are: t « 3eF, t± « Q.AeV, ti — 
0.1 — 0.3eV. ti cannot be estimated accurately and 
the limits of the range quoted are from experiments 
and band-structure calculations respectively. The crys- 
tal structure yields an onsitc potential difference between 
stacked and unstacked site of A = 0.01 eV. But as noted 
inversion through the stacked atom in each layer is pre- 
served. 



III. MEAN FIELD ANALYSIS 

We consider states that break time reversal symmetry 
through orbital loop-currents. In Appendix El we discuss 
why symmetry breaking states generated due to Vnn need 
not be considered even though Vnn ~ 2V^„„. Therefore 



Ou,j = (iclci.+H.C.y (4) 

The one particle terms can be dropped. Now a mean-field 
approximation is made with Vij{Ou.j) = r for next near- 
est neighbor bonds This leads to a time-reversal 
breaking Hamiltonian for electrons with complex hopping 
between sites: 

Uj + (Vi{clcu^j + H.C.)^ 

= Uj + inj = Uj exp^*'^''^', (5) 

for r <^t, as well as terms in energy of the form r^/21/. 
r may now be determined by minimization of energy. 

There are four kinds of states which can be generated 
through Vnnn which break time-reversal without break- 
ing translational symmetry. They can be classified by 
the signs of the fiux within different areas (triangles) en- 
closed by the next-nearest-neighbor bonds. (1) Within 
each layer, the triangles of A sublattice and B sublattice 
can have the same or opposite signs of fiux (see Fig. 
(2) Between two layers, the triangles of the unstacked 
atoms (Bl, A2) centered at the stacked atoms (A1,B2) 
can have the same or opposite signs of fiux as well [see 
Fig. |31[a-b,d-e)]. Of these the two states which break in- 
version through having flux in the stacked and unstacked 
triangle of atoms in a given layer in opposite direction 
[Fig.mb)] always have a higher ground state energy than 
the other two constructed from Fig ^a) and will not be 
considered further. We are then left with two possibil- 
ities, with flux in the stacked and unstacked triangle of 
atoms in a given layer in the same direction. They can 
be further classified by the second condition, and are 

(i) The Haldane or AHE state of the bilayer in which the 
orientations of the flux in the loops atop each other in 
the two layers are the same. The fiux patterns are shown 
in Fig. [31Ja-c). 

(ii) The Magneto-electric (ME) phase with the opposite 
orientations of the flux in the loops atop each other in 
the two layers. Fig. [31[d-f). In both cases the net fiux 
through a unit-cell is zero. 

We calculate the ground states energy in mean-field in 
the space of all the four bands. We define r of Eq. ([5]) by 
using hiitV = Vnnn and (i, j) nearest neighbors for both 
stacked and unstacked atoms in the same layer. We find 
the leading dependence in energy on the order parameter 
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(a) 




FIG. 2: The pattern of net fluxes in one layer (the top layer 
is illustrated as an example). The red and blue colors repre- 
sent the positive (counterclockwise) and negative (clockwise) 
fluxes, respectively. For each layer, there are two possible 
loop-ordered states generated through V„„„: the triangles en- 
closed by bonds connecting Al, Bl sublattices and with the 
same center have the same (a) or opposite (b) signs of fluxes, 
(a) is corresponds to the Haldane's AHE flux pattern for a 
single-layer honeycomb lattice. Here, the inversion symmetry 
is preserved because the center of the honeycomb lattice is the 
inversion point. However, the inversion symmetry is always 
broken for (b). 




(b) k (c) 





FIG. 3: The pattern of fluxes in the bilayer as viewed from 
above, (a-c) represent the AHE state while (d-f) represent 
the ME state, as discussed in the text. For each state, the 
flux patterns for the top layer (a and d), the bottom layer 
(b and e), and the net fluxes (c and f) are shown. The color 
representation is the same as in Fig (2] In (f), we use light 
red (blue) to represent weaker strenghes of fluxes compared 
to the regions with regular red (blue) colors. 



r for the AHE and the ME states to be, respectively. 



E{r)ME/t 

E{r)AHE/t 



1 



2V, 



nmi 

1 



9.056 + 228. Sr" 



9.422 + 306.4r*. (6) 



Eq. m gives the interesting result that for Vnnn ^ 0.05f , 
the magneto-electric phase has a lower ground state en- 
ergy in such mean-field calculation. Eq. IH] also suggests 
that as temperature is decreased, the first state to arise 
through a transition with Ising symmetry is the AHE 



state (with a first order transition), and then to the ME 
state at a lower temperature. We calculate the tempera- 
ture dependence of the ground state energy, from which 
we estimate the transition temperature to the ME state 
is Tc ~ 10~^t (for Vnnn — 0-U). The transition tem- 
perature to the AHE state is 5% higher. The value of 
the transition temperature is expected to be depressed 
from these estimates by fluctuations. The one-particle 
excitations in AHE state have a direct gap while in the 
ME state have an indirect gap with features at higher 
energy, shown in Fig. [S] which should be looked for in 
experiments. 



IV. NUMERICAL ANALYSIS 



To examine the property of the time-reversal symmetry 
breaking state with an applied magnetic field, we diago- 
nalize the mean-field Hamiltonian with a magnetic field 
on a finite lattice of the size of the magnetic translation. 
We calculate both the density of states and the various 
topological properties of the AHE and the ME states, 
such as the edge states and the Chern numbers, when 
they exist. The details are given in the Appendices. The 
important results are stated here. 

For the AHE state, any finite order parameter Vnnn 
will open up a gap A « 6\/3r. As verified by the edge 
modes, see Fig. |31 it has a quantized Hall conductance 
at zero magnetic field. We calculated the Chern numbers 
for the 4 bands, which are C = (0, 2, —2, 0) in sequence of 
energy. This gives a quantized Hall conductance a^y = 
4e^//i for r < 0{t±). For r > 0{t±), the Chern numbers 
are C = (2,0,0,-2) instead. 

The ME state has an indirect gap (see Fig. 0]) only for 
r > 0{t±). For small r, one obtains a semi-metal. The 
insulating ME state is not accessible within the two-band 
low-energy approximation. As shown in Fig. |4l the edge 
modes do not cross; so the ME state does not carry a 
net Hall current at zero field. The direct gaps produce 
further sharp features in one-particle excitations which 
should be looked for to confirm the existence of this state, 
as shown in Fig. [S] 

We expect that the density of states of the ME state at 
the edge of the gap becomes much sharper because the 
excitonic effects are not included in our calculation. 

We now examine the property of ME state in a finite 
magnetic field. The density of states is shown in Fig. [S] 
Unlike the AHE state, where the lowest Landau level is 
pinned to the edge of the gap [2l|, the lowest Landau 
level(s) shifts up in energy, leading to a larger gap. The 
magnetic field dependence of the gap is shown in Fig. |6l 
and agree with experiments when extrapolated to low 
fields of the experiments. 
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FIG. 4: The bandstructures with edge modes (assuming an 
open boundary condition along the zig-zag edge) for AHE 
(top panel) and ME (bottom panel) states. For the non- 
interacting bilayer graphene, there are four bands: two low- 
energy valence and conduction bands touch each other at two 
momentum points K and K'. The other two high-energy bands 
are split with a scale of t^. In the AHE state, band gaps 
are opened at both K and K' due to the order parameter r, 
A ~ &\/ir. In the ME state, there are also direct band gaps 
at K and K', ~ 3\/3r. But the center of the direct gaps 
shift in different directions of energy for K and K'. When 
Ad i±, an indirect gap opens. For illustration purposes, we 
take = 0.2t and r = 0.04t here, which are slightly bigger 
than realistic values. 



V. CONCLUSION 

The principle experimental observation of an insulat- 
ing state at charge neutrality in bilayer graphene, with 
a gap growing monotonically with a magnetic field, is 
the signature of a novel state of matter. We have shown 
that these results are only consistent with a state that 
breaks both the time-reversal and inversion symmetries, 
but not translational invariance. When a Haldane phase 
is realized in each layer, the even parity combination of 
two layers (i.e. AHE state) has edge modes and finite 
conductivity, while the odd parity combination (i.e ME 
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FIG. 5: The density of state of ME state without and with a 
magnetic field. Here t± = Q.2t and r = 0.04t. The zero-field 
DOS has been multiplied by a factor 5. The magnetic field 
strength is specified by the fiux per unit cell <j) = 2n/M where 
M is an integer (see Appendix B for more details). 
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FIG. 6: The gap as a function of the external magnetic field. 
Here t± — 0.2t. The magnetic field strength is in terms of 
the flux per unit cell (j) ~ 2n/AI, and M has been chosen as 
128,256, ...,2048. 



state) does not. We have shown that the ME state is 
the ground state of interacting bilayer graphene, and its 
properties are in agreement with those of the observed 
insulating state. We suggest conductance experiments at 
higher biases to confirm the ME state through observa- 
tions of the jumps in the density of states of Fig. [Sj 

We have also considered the possibility of the actual 
realization of the Haldane AHE state in a single layer 
graphene by our mean-field procedure. Again the nearest 
neighbor interactions are ineffective and need a very large 
value to get any change in symmetry but a next nearest 
neighbor interaction Vnun approx larger than about O.li 
leads to the Haldane state. This requires a dielectric 
constant to be smaller than about 10, which appears to 
be close to what experiments report. 



5 



Acknowledgements 

We acknowledge extensive discussions with Professor 
C.N. Lau. L.Z. also acknowledges discussions with Pro- 
fessor Donna Sheng. This research was partially sup- 
ported by NSF grant DMR-0906530 (L.Z. and C. M. V.), 
and UCR initial complenient(V.A.). 



Appendix A: Ordered Loop Current States 
generated through Vnn 

For the nearest neighbor interaction Vnn , the generated 
loop currents flow around the perimeter of a hexagonal 
cell with net flux. If we allow flux (positive or negative) 
in any unit-cell, we run into the problem of frustration 
of an Ising model in a triangular lattice. This in gen- 
eral is not a possible ground state because unlike with 
spins, alternate arrangements with lower energy are pos- 
sible. The simplest is a ^/3a x -y/Sa structure of hexagons 
depicted in Fig. [7] in which a hexagonal cell with zero 
net flux is surrounded by six hexagonal cells, each with a 
Kekule pattern of currents in the links with alternate net 
positive and negative flux. This state breaks the transla- 
tional symmetry while the unit cell is enlarged to include 
3x3 unit cells of a translational invariant honeycomb 
lattice. 




FIG. 7: A possible ordered loop-current state generated 
through Vnn- Only one layer is shown. The red, blue and 
white colors represent the net positive (clockwise), negative 
(anticlockwise), and zero fluxes in a hexagonal cell. 

We assume that the flux pattern in the top layer follows 
Fig. [3 which takes a sequence (+,-,0) from left to right. 
The bottom layer could follow the same pattern, or takes 
another pattern (-,-|-,0). Therefore, there are two types 
of ordered loop current states (the other combinations 
are equivalent to either of these two states by rotation). 
We carry out a mean-field calculations on the ground 
states energy for these states, and find that the state 
with different patterns for two layers are lower in energy. 
This lower energy state is allowed only if Vnn ^ 2.06i, 
i.e., be comparable to the bandwidth. The reason is that 
such a state does not have a "nesting" periodicity and 
therefore does not use uses any singularity in the joint 



density of states. 



Appendix B: Numerical Methods 

The mean-field Hamiltonian of the ordered loop cur- 
rent states can be readily diagonalized for each momen- 
tum point (fca;, ky). However, to examine their properties 
under an applied magnetic field, as well as the topologi- 
cal properties such as the edge modes, we also carry out 
a real-space calculation. 




FIG. 8: The bandstructures with edge modes (assuming an 
open boundary condition along the zig-zag edge) for AHE (top 
panel) and ME (bottom panel) states. Here we take t± = 0.2t 
and r = O.Oit. 



Following a common numerical practice for honeycomb 
lattices, we describe the bilayer graphene by a Nx x Ny 
lattice, with four atoms on each lattice site (see Fig. [S|). 
The real-space coordinations of each atom (and the mo- 
mentum) in the orthogonal x-y axis can be easily ob- 
tained from this "deformed" -lattice representation. 

The effect of an external magnetic field can be captured 
by imposing a phase to the hopping term to each bond 
Uj tijc^^ We adopt the Landau gauge Ay = 
~Bx to take advantage of the translational invariance 
along y-dircction. For the periodic boundary condition, 
we commonly take the flux due to the external magnetic 
field in a unit cell to be = 27r/M, where M is an 
integer proportional to N^- M — 1000 is equivalent to a 
magnetic field strength B w 30r for graphene systems. 

In general, we choose periodic boundary conditions 
along both directions ( a torus). This provides a verifi- 
cation of our momentum-space calculations for the zero- 
field. For the edge modes calculation, we choose an open 
boundary along x-axis, which corresponds to a cylinder 
with the zigzag edge. 
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Appendix C: Additional Results 
1. Density of States 

In Fig. m we show the density of states (DOS) for the 
ME state with various values of the order parameter r. 
We also show the result for the non-interacting case as 
well. For r = 0, the DOS at w = is finite, which is due 
to the parabolic dispersion E ^ k^. There is a sudden 
jump of DIS at a; « zLt±/2, which is the energy scale of 
the gap between top two valence bands (bottom conduc- 
tion bands). When r is small, the conduction band at K 
point and the valence band at K' point have overlaps in 
energy. Therefore, it remains a semi-metal state. A full 
gapped state happens when r > t^/(3\/3). 




I ' ■ ■ ' ' ' ' 

-0.4-0.3-0.2-0.1 0.1 0.2 0.3 0.4 



CO 

FIG. 9: The density of states for the ME state. Various 
values of the loop current order parameter are shown. The 
non-interacting case r = is also shown in comparison. Here 
t± = 0.2t. 



2. Chern numbers 

In the main paper, we have shown the topological prop- 
erties of the AHE and ME states by showing their edge 



modes, calculated by assuming a periodic boundary con- 
dition along y-axis and an open boundary condition along 
X-axis. Judging from the crossed edge modes, we learn 
the AHE state indeed has a finite Hall conductance at 
zero-field. Similarly, the ME state does not have any 
Hall effect at zero-field. 

This topological property can also be checked from the 
Chern numbers. The Chern number for a given band is 
defined as 

f <fk ^ a^*(k) ^ 9^n(k) 

where ^'„(k) is the wavefunction for n-band. The Hall 
conductance is given by (Txy = Ce^/h. We calculate the 
Chern numbers in a momentum-space algorithm, by sum- 
ming the Berry curvature in small areas {Ak^, Aky). 

For the AHE state, any finite order parameter r will 
open up a gap. We find the Chern numbers for the 4 
bands are C = (0, 2, —2, 0) in sequence of energy, i.e, the 
valence and conduction bands away from fermi energy do 
not carry any Chern numbers. However, when r > t^, 
the Chern numbers become C — (2,0,0,-2) instead. In 
either case, a finite Hall conductance is a robust feature. 

For the ME state, an indirect gap is open only for a 
finite order parameter r > t_\_. However, as long as the 
four bands do not touch with each other (having direct 
band gaps), the Chern number calculation is robust. We 
find the Chern numbers for the four bands of the ME 
state are all zero, indicating a none-quantum Hall state. 
When r S> we find the Chern numbers become C = 
(1,-1,1,-1). In this limit, each layer is a Haldane state 
while t± only acts as a perturbation, slightly splitting the 
stats from two layers. Still, there is no Hall effect near 
the Fermi energy. 



3. Edge Modes in a finite magnetic field 

In the main paper, we have shown the edge modes for 
AHE and ME states at zero-field. For completeness, we 
also show the edge modes in a finite external magnetic 
field, in Fig. [TUl 
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